General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



AgRISTARS 


“Made available ynder NASA sponsorslilp 
in the interest of early and wide dis- 
semination of Earth Resources Survey 
Program information and without liabUHy 
for any use made thareot.” 

Conservation and Pollution 


E83'1042i 

CP-53-04417 _ 

A Joint Program for 
Agriculture and 
Resources Inventory 
Surveys Through 
Aerospace 
Remote Sensing 


THE CONTINUOUS SIMILARITY MODEL OF BULK SOIL-WATER EVAPORATION 


(E83- 10426) THE COHTIIIOODS SIHILIfilTI BODEi 
OF BULK SOIL-IATBB EfAPOBATIOil (Haryland 
Univ.) 55 p ac A04/BF A01 CSCL 02C 


G3/43 


B83-34411 

Unclas 

00426 


ROGER B. CLAPP 


REMOTE SENSING SYSTEMS LABORATORY 
DEPARTMENT OF CIVIL ENGINEERING 
UNIVERSITY OF MARYLAND 
COLLEGE PARK^ MARYLAND 20742 




LynNon B. Johnaon SpeoeCanlef 

Houston, Texas 77058 


original par® 

^llAl 1 1 ■ 


TCCHHICAt RIPOftT STANDARD TITLf PACE 


1. Ra^ot Na 




2. C***ra<*«M AtctttiM 


The Continuous Similarity Model of Bulk 
Soil-Water Evaporation 


3. R»ci^i«n«‘t C*mI *9 N*. 


$. R»p«r« Dal* 

April 1983 


6. Rarbrminf Orgmitatien Cadt 


7. A»tlior(») 


Clapp, Roger B. 


B. Rarlerming Organitofian Rapart Na. 


9. Railarming Organ! laliaa Nama and Addratt 

Remote Sensing Systems Laboratory 
Department of Civil Engineering 
University of Maryland 

Park. MD 2Q742.. 


12. S^ent«ri«i 9 Af#ncy Nom# ond Addftti 

USDA-ARS 

Beltsville, MD 20705 


15- Notes 


^16- Abitroct 

I See attached. 


10. Wark Unit Na. 


11. Caniroct ar Grant Na. 




13. Tjrpa a! Rapart and Pariod Covarad 


14. Spantaring Agancy Cada 


17. Kay Wardr (S- 'acted by Aulhaf(a)) 

Soil-water dynamics, evaporation, 
hydrological modeling. 

18. Oittfibatian Statamant 

UNLIMITED 

y-5 

19. Sacarity Clattif. (of rhit rapart) 

20. Sacarity Clai alf. (of thit page) 

21. Na. af Pagai 

22. Price* 

Unclassified 

Unclassifed 

50 



For nalc by the Clrarin|thouar for Federal Scieniiftc and Technical Infotiaalion, Sgrin|tfield, Virginia 22I5I- 











THE CONTINUOUS SIMILARITY MODEL 
OF BULK SOIL-WATER EVAPORATION 


BY 


Roger B. Clapp 



ABSTRACT 


In the past, evaporation from unvegetated soils has been sinulated by 
using complex, numerical models of simultaneous moisture/heat flow and by 
using event models. The latter are attractive due to their relative simplicity 
although they have been difficult to apply to field conditions. This is so 
because variable climatic conditions and redistribution of soil water following 
infiltration both complicate the evaporation process. 

Both factors are incorporated into the continuous similarity model of 
evaporation. In it, evaporation is conceptualized as a two-stage process. 

For an initially moist soil, evaporation is first climate-limited, but later 
it becomes soil-limited. During the latter stage, the evaporation rate is 
termed '*evaporability," and mathematically it is inversely proportional to 
the evaporation deficit. In this relationship the model is similar to the 
Green-Ampt model of infiltration in which the infiltration rate is inversely 
proportional to cumulative infiltration. Unlike the Green-Ampt model, the 
evaporation model requires numerical integration. It also includes a functional 
approximation of the moisture distribution within the soil column. 

The model was tested using data from 4 experiments conducted by the USDA-ARS 
near Phoenix, Arizona; and there was excellent agreement between the simulated 
and observed evaporation. The model also predicted the time of transition to 
the soil-limited stage reasonably well. 

For one of the experiments, a third stage of evaporation, when vapor 
diffusion predominates, was observed. The occurrence of this stage was related 


to the decrease in noisture at the surface of the soil. The continuous 
sinilarity nodel does not account for vapor flow. 

The results of this study show that climate, through the potential 
evaporation rate, has a strong influence on the time of transition to the 
soil-limited stage. After this transition, however, bulk evaporation is 
independent of climate until the effects of vapor flow within the soil 
predominate. 
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Definition sketch. For the CSM, E* corresponds to the area 11 
within the bold line. In time, E* in increased by evapora- 
tion (hatched area) and decreased by redistribution (gray 
area). The triangles indicate the depth of drying. Because 
the changes in E* due to evaporation and redistribution are 
evaluated continuously in the model, the overlap between the 
incremental changes shown in the diagram presents no problem. 


This procedure calculates the rate equations of the state 16 
variables; it is incorporated within the numerical integra- 
tion scheme. 

Calculated and measured evaporation. The circled symbols 24 
Indicate data used to generate the average PE rate. 

Calculated and observed moisture profiles for the March 27 

experiment . 

Model variables. The bold lines indicate E* and on 29 

the appropriate axes. Symbols indicate the observed E* 


for March, and the hatched areas indicate the depth interval 
of maximum observed 6 where the actual depth of drying 
probably occurred. 

Transformed evaporation data. Each data symbol corresponds 34 
to one day, and certain symbols are individually identified 
by the number of days following irrigation. Simulated 
results for July, September, March, and December extend to 
7, 14, 37 and 14 days, respectively. Circled symbols 
indicate data used to generate the average PE rates. 

The daily average moisture content based on measurements 36 
from the uppermost 5-mm layer. 
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PREFACE 


This report is the sequel to a previous AgRISTARS publication entitled 
’’The Desorptivity Model of Bulk Soil-Water Evaporation.’* As with the previous 
one« this report is not concerned with the nechanics of remote sensing, rather 
it deals with the physics of the near-surface region of the soil. 

Of all the relationships developed herein, the item that may be most 
immediately helpful to the problem of interpreting moisture measurements 
obtainable by remote sensing is the functional approximation to the soil 
moisture profile described in Section 2. This approximation applies only 
when the soil surface is relatively dry, and it estimates the 6 - z relation- 
ship only down to the zero-flux depth. Since this depth varies in time and 
rarely is expected to be greater than 50 cm, the approximation itself 
characterizes the difficulty in extrapolating the moisture distribution from 
surface measurements. Recognition of this difficulty reinforces the veiw point 
that simpler models — like the one offered herein, and as opposed to cooqjlex, 
numerical ones — offer the best avenue for infusing both soil physics and 
remote sensing into hydrological modeling. 


Roger B. Clapp 
March 1983 
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CHAPTER 1 


INTRODUCTION 

Physically based aodels of avaporation frcHs unvogetated soils can be 
divided into essentially two grotqps: nunerical node Is of sinulteneous 

aoisture/heat flow (e.g., Hilly, 1982) and event nodels (e.g., Gardner and 
Hillel, 1962; Gardner et al., 1970 Stb). Models of the first group lire con^lex 
and data*intensive, hence they are expensive to r\m and sometimes they give 
results that ere difficult to interpret. The sheer volioe of such results 
can possibly obscure sisg>le nechanisas that cm be controlling the evapora- 
tive process. In contrast, event aodels of evaporation are attractive due to 
their simplicity although their ^plicability is liaited by their tmderlying 
assuaptions and by the fact that they have never been satisfactorily field 
tested. Given these differences between aodels, it would be desirable to know 
what aspects of evaporation can be explained by a simple event model. From 
a practical stan^oint, it is desirable to have a field-tested event aodel that 
could be aerged with an infiltration model to yield a simplified, yet physically 
based model of the field water cycle. In turn, such a aodel could be applied 
to large-scale problras where the inherent uncertainties favor a siiq>lified 
modeling mpproach. 

Event models of evaporation simlate bulk flow, i.e., the rate of evaporation 
averaged over 24 hours; and they assume that evaporation can be divided into 
distinct stages. For an initially moist soil, the evaporation rate is mostly 
controlled by meteorological conditions. In time, a transition occurs when the 
rate becomes limited by the soil's ability to transport water to the near-surface 
region where evmporation actually occurs. The two intervals are identified as 
the climste-liaited and tte soil-liaited stages, respectively, or simply at 
stages 1 and II. 
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As tht soil bee<»«s very dry there is e possibility of a third stage 
When evaporation is aaintained at a low« and perhaps steady, rate by the dif- 
fusion of water vapor within the soil. Both the occurrence and the iaportance 
of stage- 1 II evaproation tinder field conditions are open to question. In 
addition to these stages, Idso, et al. (1979) described a transitional interval 
between stages I and II tdien the soil surface exhibits both wet and dry areas. 

This "patichiness'* is probably due to spatial variability in soil properties 
and/or Bicroclinate. The uodel developed herein does not account for this varia- 
bility or its effects on ev^oration, although this phenoeenon was observed 
in the experiaents that provide the data used to test the sodel developed herein. 

The fundaaental assuaption of event aodels is that stage- I I evaporation 
is essentially independent of weather conditions. Nevertheless, it is coamon 
to see models that link stage-ll evaporation to meteorological variables. For 
instance. Barton (1979) ad^ited the equation of Priestley and Taylor (1972) for 
potential evaporation to account for water loss from unsaturated surfaces: and 
Beese, et al. (1977) related the actual evaporation rate to the potential rate 
using an eipriicism based on measures of soil-water- suction. Clearly, there is 
a basic contradiction between these approaches and the event models. For 
periods less than 24 hours, meteorological conditions are i^>ortant in controlling 
evaporatiw) because instantaneous rates of evaporation are correlated with net 
radiation, with the highest rates often occurring near solar noon. However, the 
focus of this research is on bulk evaporation, defined as the flux to the 
atmosphere averged ever 24 hours. Ihe question is whether bulk evaporation 
under field conditions can be adequately described by an event model. 


3 


Thii qiMStion was afflnttivtly in tht prtvious paptr (CUpp, 

1983, hmctforth rtferrtd to as PAXt I) in which stAfe>II tvi^orAtion was 
dAScribod by the dtosorptivity aodel. Despite its Accurecy, thAt aof*el hAS two 
sifnificAnt drAwbAcks: (me prActicAl, the other conceptuAl. First, it requires 

An if^ependent estlMte of the tisie of trensition to stAge II. Although this 
tistt CAn be estisAted a po»tsriori from chsnges in Albedo, no predictive 
nethod WAS given. Second, the asin equetion of the SK>del requires s fixed value 
for desoxptivity, a key perAseter in the aodel that is dependent on, aaong other 
things, t)w Bioisture content within the soil coluan. Due to the redistribution 
of Boisture that follows infiltration, the noisture Available for evaporation 
decreases in tiaw so that the desoxptivity also decreases. Consequently, the 
evaporation sK>del includes an esq^irically chosen nethod for generating a 
r^resentative value for desoxptivity. The substitution of a physically based 
relationship for this eapirical averaging technique would constitute an isprovesient. 

1.1 Objectives 

With these probleas in aind, the first and second objectives of this paper 
are to describe and to test a aodel of bulk evaporation that estiaates the tiae 
of transitim and that incorporates the effect of redistribution directly. The 
new aodel, teraed the Continuous Siailarity Model (C5M), was tested using the 
saae data used to develop the desoxptivity aodel in Part I. The success of the 
aodel test provides aore evidence that bulk evaporation during stage II is 
essentially on isotheraal process. In other word*, idten evaporation is soil- 
Halted the effects of tesperature and energy on the rate of evaporation -- when 
averaged over 24 hours — are negligible. 

Hie aodel test also provides indirect evidence that tiM assuaptiim of 
distinct stages is valid; but the ability to aatch cuaulative ev^oration over 
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• period of, eay, 14 days does not prove ^he existence of separate si jes. 
Consequently, the third objective it to exanine the evaporation 4«ta directly 
to assets hw distinct «id identifiable the stages of evaporation are. Of 
prise isportance is the occurrence of stage III when vapor flow px^dmainates . 
Because the CSM does not account specifically for vapor diffusion, sodel errors 
in the calculated evaporation night occur under very dry cemditions. 

The fourth objective is to relate the observed stages of evaporation to 
the aeasured noisture content at the surface of the soil. Although the CSM 
does not require the surface s^isture as an input, the surface noisture is 
is^ortant for a variety of reasons. For instance, it is strongly related to 
albedo and thus to the atmospheric energy budget. In addition, surface mois- 
ture can be measured indirectly by remote sensing; and it would be highly 
desirable to use such measurements to extrapolate the smisture profile, i.e., 
the moisture distribution with depth. Among other things, the includes 
a itmctional approximation for the moisture profile (appropriate when the 
surface is relatively dry), and thi^ 'Sion may aid in the interpretation 
of reMtely sensed data. It follows that the sc(^e of this report ranges from 
the specific, conceptual development of the continuous similarity model to a 
broader investigatim of the process of soil-water evaporation. 



CHAPTER 2 


MODEL DEVELOPMENT 


Tlte three sain steps in the development of the continuous similarity model 
are presented in the next three sections. In the first section, the theoretical 
basis of the model, which was described in detail in Part I, is briefly reviewed. 
The fundamental expression for the rate of evaporation during stage II is intro> 
duced and then rearranged so that the rate is proportional to the reciprocal 
of the evaporation deficit, E*. This inverse relationship was first explored 
by Philip (1957), and in this form the evaporation model is analogous to the 
infiltration model of Green and Anpt (1911) in which the infiltration rate is 
inversely proportional to the cumulative infiltration. However, as is shown 
in the second section, E* in the is not singjly the cumulative evaporation. 
Although E* increases as evaporation continues, it is modified by the effects 
of redistribution. Quantification of the redistribution effect requires a 
functional approximation to the moisture profile, and this approximation is 
developed in the third section. Whereas the Green-Ampt model assinnes a 
simple rectangular shape for the moisture profile, during evaporation the mois- 
ture profile is expected to be curved; this nonlinearity is critical to the 
estimation of the evaporation rate. 

In the fully developed model, evaporation is described by the behavior 
of two dynamic state variables: E* and 6^, the moisture content at depth. 

The temporal change in 6^ results from the redistribution process which 
under certain circianstances can be assumed to be independent of evaporation. 
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In addition, the fully developed aodel has two special attributes. First, 
unlike the Green-An^t model, the CSM requires numerical inteRration to advance 
the solution throuRh tisw. Second, the has two forms correspondins to 
the two commonly used expressions for soil -water diffusivity. 

2.1 Desorption and the Evaporation Deficit 

Part 1 describes the desorption problems of soil water dynamics and the 
associated "similarity solution." Together, these mathematical relationships 
provide the theoretical basis for the CSM. The desorption problem represents 
evaporation from a semi-infinite column of homogeneous soil having a uniform 
initial moisture content. The surface of the soil is assumed to be dried 
instantaneously and thereafter maintained at a low, steady moisture content. 

In addition, drainage effects are assumed to be negligible. At all times, the 
evaporation rate is soil-limited. In Part I, it was shown that the rate is 
sensitive to the soil's diffusivity funcr.ion, as well as to the initial mois- 
ture content. However, it was also shown that the evaporation rate is largely 
indpendent of the surface moisture content provided that the moisture content 
is less than a critical value. 


For the desorption problem the similarity solution is essentially a 
mathematical description of the moisture profile. For any diffusivity function, 
there is a unique function: 

8(y) or 0 ( 2 , t) 

describing the moisture profile, where 6 is the volianetric moisture content 
and y is the similarity variable defined by the relationship: 


i(t - V 


-i/j 
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The parameter is an en;>irical time delay that accoimts for the initial, 
climate- limited stage of evaporation; it is not included in the rigorous 
solution of the desorption problem. As shown in Part I, integration of the 
similarity moisture profile with respect to z yields the evaporation rate: 

dEj/dt « 1/2 A(t - (1) 

which is the fvindamental expression in the desorptivity model. In Eq. (1) 

E 2 refers to evaporation during stage II, and k is the desorptivity parameter 
which is functionally dependent on the moisture content at depth 6^. The 
latter parameter is equivalent to the initial moisture content in the desorption 
problem; and in Part I, it was shown that 6^ and therefore A decrease in time 
after the soil is initially moistened during infiltration. The change in 6j 
is caused by redistribution, and as mentioned before, this change has a signi- 
ficant effect on the evaporation rate. 


In the past, analytical solutions for 6(y) and A were unavailable for 
commonly used diffusivity functions, but in Part I two approximations for 
A were presented. They were developed for two alternative expressions 
for diffusivity, specifically: 


and 


D,(e/ep 


exp (ae) 

referred to as power D and exponential D, respectively. The parameters D , D , 

9 O 

c and a are usually fitted to D-e data by regressioiv and there is no substantial 
evidence as to which expression for diffusivity is more appropriate. The 
approximations for A are 
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[ 12 D 


e ^ 
s s 

I)(c*4) 


e, c*2l 


1/2 




C2) 


for power D, and 


A 

e 


11.3 expCaSp 
<nr(a0. ♦ 1.85) 

_ ^ m 


1/2 


(3) 


for e 3 q>onential D. Both expressions are independent of the surface moisture 
content, 6^, and they are considered to be accurate where 0^ < 0^, the 
critical moisture content (as defined in Part I). 


In the new evaporation model, Eq. (1) for the evaporation rate is modified 
so that the independent variable t is eliminated, and the rate is calculated 
from the cumulative evaporation — in much the same way as the infiltration 
rate in the Green-Ampt equation is calculated from the cumulative infiltration. 
To eliminate t, first (1) is integrated with respect to (t - t^): 

E* « ^(t - (4) 

where E* is the evaporation deficit. At this point, E* is equivalent to 
the cumulative evaporation; however, a more general definition is given in 
the next section. Next, (4) is rearranged and substituted back into (1) to 
yield the expression: 

dE2/dt - U^/2)/E* (5) 

which is merely a restatement of (1). 


Under field conditions, the similarity solution cannot be applied 
rigorously because simultaneous redistribution decreases 0^ through time. 
Despite this redistribuiton effect, the key supposition of the new model is 
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that (5) and a corresponding similarity noistnre profile, 6(y), are contin- 
uously valid through stage II, hence the name ^'continuous similarity model." 
For simultaneous redistribution (5) is more appropriate than (1) because 
(5) relates dE^/dt to the instantaneous distribution of the moisture within 
the soil column when the value of A in (5) is the instantaneous value. On 
the other hand, A in (1) must be fixed, so when redistribution changes 6^ 
a tin»- averaged A must be specified. There is no physical basis for any 
particular averaging scheme. 


For simplicity, (5) can be rewrtten as 

dE2/dt - «Cej)/E* (6) 

2 

where ^ ^ A 12 and e is functionally dependent on 6^. From (2) and (3) it 
follows that the approximations for « are: 

6 D ef 

e . s s 

P (c+l)(c+2)ir 

for power D, and 

5.65D e, exp (a6,) 

♦, 2 .^ — ( 8 ) 

air(aej ♦ i.ss) 



for exponential D. • represents the capability of a homogenous soil to 
move liquid water to a dry surface. Because (7) and (8) are singly manipu- 
lations of the approximations for A, the critical moisture content still serves 
as a guideline indicating how dry the surface must be. Again, for 6^ < 6^ 
the evaporation rate determined from (6) using (7) or (8) is considered to 
be a good approximation to complex, nximerical approximations of the original 
similarity solution. The effects of the surface moisture on the evaporation 
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rate are nore fully ej^lored in the discussion section. It is sufficient 
at this point to note that 6^ values tend to be quite large* certainly 
larger than *'air-dry*' values of 0. 


2.2 Effects of Redistribution and Climate 

Following infiltration, the wetting front continues downward due to 
gradients in both the natric suction and the gravitational potential. Because 
no moisture is supplied to the surface, the moisture content behind the front 
decreases. The rate of decrease, termed the redistribution rate is initially 
large but decreases through time. I hypothesize that E* is conditioned by 
redistribution, i.e., E* simultaneously increases due to evaporation loss at 
the surface and decreases due to redistribution, as illustrated by the defi- 
nition sketch in Figure 1. From that figure, it can be seen that 


I. 

E* « J °(e - e)dz 
0 ^ 

and that the rate of change of E* is given by the expression: 


(9) 


dEVdt - dE/dt ♦ z^Cde^/dt) (10) 

where dE/dt is the actual evaporation rate, and is termed the depth of 
drying. Technically, z^ is the depth to the zero-flux point within the soil 
column; but for simplicity, in the CSM z^ and 6^ are specified at the point 
of maximum e above the wetting front. It follows that z^ increases in time 
as evaporation removes water from progressively lower depths. The rate of 
redistribution, d6^/dt, is inherently negative; thus the redistribution 
effect tends to decrease E*. From (10) it can be seen that if de^/dt ■ 0, 
as required for the rigorous similarity solution, then E* is equivalent to 
the cumulative evaporation. 
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FIGURE 1. Definition sketch. For the CSM, E* corresponds to the area 

within the bold line. In tine, E* is increased by evaporation 
(hatched area) and decreased by redistribution (gray area) . 

The triangles indicate the depth of drying. Because the 
changes in E* due to evaporation and redistribution are evaluated 
continuously in the model, the overlap between the incremental 
changes shown in the diagram presents no problem. 
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To simulate evaporation as a two-stage process, the time of transition 
to stage II is computed isplicitly. This is done by calculating dE^/dt from 
(5) and by applying the following relationships: 

dE2/dt ^ PE dE/dt « PE, 

dE^/dt < PE dE/dt ■ dEj/dt. 

These relationships are analogous to those of Mein and Larson (1973) for 
determining the change from climate-limited to soil-limited infiltration. 

Hillel (1971) called the soil-limited rate of infiltration the soil’s 
"infiltrability;” thus for the sake of symmetry, dEj/dt is termed the soil's 
"evaporability.” Evaporability serves to indicate both the beginning of 
stage II and the evaporation rate during stage II. Evaporability depends 
on E*, and E* is calculated from, among other things, the depth of drying. 

The estimation of requires a mathematical expression for the moisture 
profile, and that is developed in the following section. 

2.3 Moisture Profile 

The specification of the moisture profile requires an alternative 
approximation to the similarity moisture profile. The alternative approach 
is based on the Kirchhoff transformation which was first applied to the 
process of soil-water evaporation by Gardner (1959). The transformation serves 
to siiq>lify the governing equations for the desorption problem, and the reader 
is advised to compare the following equations with those of the original 
desorption problem in Part I. With the Kirchhoff transformation defined as 

0 0 

U - D de/a^ ^ D de (11) 

®0 ®0 

and substituted Into Eq. (1) of Part I, the partial differential equation 
for soil-water dynamics (without gravity flow) becomes 
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It ■ 


( 12 ) 


With this transformation the initial and boundary conditions of the desorption 
problem becone 


t«0, z>0, U«1 
t>0, z*0, U»0 
t > 0, z U ■ 1. 

From these transformed equations it is possible to derive simple functions 
for the moisture profile applicable whenever the surface is relatively dry. 


Eq. (12) has an approximate solution; 

U ■ erf[z/(4 D** t)^] (13) 

where erf is the error function and D** is an average diffusivity dependent 
on U. (Note that D** is not readily related to the mean weighted diffusivity, 
D*, of Part I.) In this solution, 6 approaches 6^ asymptotically so (13) 
does not computer a finite depth of drying. However, with two simplifying 
assunptions, (13) yields a useful approximation to both the moisture profile 
and Zj. 


First, the error function is approximated by a simple linear function 

erf(X) * X 

over the interval 0 < X < 1. The reason for this approximation can be inferred 
from Figure 1 in Gardner (1959) which shows that the nonlinearity in the 
confuted distribution is derived mostly from the transformation of the 
l.h.s. of (13) and not from the erf on the r.h.s. 
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Second, Eq. (12) for U can be simplified by assuming that * 0 and 
6^ ■ for power D and exponential D, respectively. These substitutions 
■ake no appreciable difference to the calculation of bulk evaporation as 
shown in Part I. Hence, for power D, substitution into (13) yields 

(O/ep*^*^ • 2/(4 D** t)** 

Because z * at 6 > 6^ the expression for the moisture profile is simply 

(e/Bj)*'''^ - 2/z^ (14) 

The development of the moisture profile through time is reflected in 
Zj which is time-dependent. However, in the model z^ is not related directly 
to time but rather to E* which also increases in time. First, (14) is 
rearranged so that 

Next, integration with respect to z from the surface to z^ and substitution 
of E* from (9) yields 

z^ - (c+2)EV6j^ (15) 

For exponential D, the identical procedure yields 

exp[a(6 - 6j)] « z/Zj (16) 

and 

Zj • oE*. 

Through substitution into (6), the governing equation for the model is either 

dEVdt - dE/dt ♦ [(c+2) E* d6^/dt]/B^ 


(17) 


15 


or 
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dE*/dt ■ dE/dt ♦ a E* de^/dt 


depending on the diffusivlty function. 


(18) 


2.4 Model Overview 

Either (17) or (18) serves as the rate equation for the state variable 
E*. Another equation is needed for the other state variable 6^, but in this 
study a single, en^irical expression is used. In the nodel test presented 
later, 6^ is described by the power function: 

* -ft 
« 0 t 

and values for the parameters 0 and 0 are given in Part I. The form of the 
redistribution function is not critical, and Gardner, et al. (1970a) derived 
several functions from physical principals. It is assumed that redistribution 
is independent of evaporation, which is reasonable where antecedent infiltra- 
tion is large; but the interaction between processes undoubtedly increases as 
antecedent infiltration decreases. This assumed independence is recognized as 
a possible source of error deserving further investigation. 

Both (17) and (18) for dE*/dt are complex; therefore E* must be determined 
by numerical integration. Although the system is uniquely described by the 
tesporal variation in E* and 0^, it is useful to integrate dE/dt to yield the 
cumulative evaporation, E. To generate the results reported in the next section, 
the rate equations for dE*/dt and dE/dt were integrated using a second-order, 
Runge-Kutta routine with a 2-hr tismstep. for exponential D, the required 
calculations to determine the rates are diagrammed in Figure 2. In the first 




FIGURE 2. This procedure calculates the rate equatims of the state 

variables; it is incorporated withJn the nunerical integration 
scheme . 


I 
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few hours of the sinulstion, dE*/dt is negstive end Urge in negnitude because 
the power function for 0^ and its derivative are not valid. During this tine. 
E* is constrained to equal zero in the aodel. From a physical viewpoint, this 
constraint is reasonable because initially no ^^reciable near-surface deficit 
is fomed. Instead the evaporative losses are distributed nearly uniformly 
with depth. 

2.5 Cca^arison to the Infiltration Model 

The model is now fully described, and the reader may choose to skip 
directly to the model test. However, it is enlightening to compare the basic 
equation for evaporability to the equation for infiltrability developed by 
Green and As^t (1911). To do this, it is first necessary to describe their 
model. 

For infiltration into a homogeneous soil having a uniform initial moisture 
content. 0^. the infiltrability is given by 

dl/dt - Kj(0j - 0p ♦ Kj (19) 

where I is the cinulative infiltration. K is the hydraulic conductivity at 

9 

saturation, and is the wetting front suction. The water pressure at the 
surface is assumed to be zero. Neman (1976) and others have derived (19) 


directly from Darcy's Law and in so doing, have defined 0^. 
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whtre k is the relttlve conductivity which is the ratio of unsaturated to 
saturated conductivities (K/K^), and is suction which is dependent on 
6 via the soil's characteristic curve. Because the integral is insensitive 
to i|)^ is effectively constant for a given soil. 

To truisfons the expression for evaporaoility into a lorn similar to 
(19) it is necessary to substitute expressions for ^ and K for diffusivity. 
This can be done using the relationships of Cai^bell (1974) for the moisture 
characteristic: 

♦ ■ ^ 

and for the conductivity function: 

where ^ is the hypothetical at 6 . By allying the conventional trans* 
formation for D, Caapbell's parameters can be related to the parameters in the 
expression for power D (as was done in Part 1). Similar relationships for 
exponential D are nut available. 

U'itg algebra and calculus it can be shown that 

dEj/dt • W <'£(0 i)/E* (HO) 

where W is a weighting factor, and is now an "unsaturated wetting front 
suction" which is dependent on via the expressioii 

tKBj) 

♦f ■ 0 ^ . 

The integral is identical in form to that describing the saturated 



19 

A tera*by>term cnoparison of expressions for evaporebility and infiltra- 
bility indicates some of the inherent differences between the two processes. 
First, it is obvimis that the constant flux term reflecting gravity flow during 
infiltration is dropped in the equation for evaporability (although drainage 
effects are calculated indirectly in CSM via the time -dependency of 6^) . 

As for hydraulic conductivity, evaporability is related to just as 
infiltrability is. At first glance, the inclusion of is counterintuitive 
because the soil system is totally unsaturated during evaporation. This 
apparent inconsistency is resolved by noting that tj»^ is inversely proportional 
to so that this factor is effectively canceled. 

The expressions for both infiltrability and evaporability contain terms 
for the moisture deficit. This deficit indicates the pore space that is filled 
(infiltration) or depleted (fviporation) during the process. For infiltrability 
the term is (6^ - 6^) whereas for evaporability it is simply 9^. In essence, 

6^ is an effective moisture deficit during evaporation because (20) applies 
even when the surface is not cosppletely dry (i.e., when 6^ > 0). 

The insensitivity to surface conditions is more evident if (20) is 
rearranged so that the evaporation rate is inversely proportional to the depth 
of drying, instead of the evaporation deficit. Substitution of (15) into (20) 
yields 

dE 2 /dt - (c+2) W «'£(ej)/ij]. (21) 

In this form, the moisture deficit has disappeared, and none of the remaining 
terms are dependent on the surface moisture. 
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Eq. (21) is analogous to the sii^list form of the Green-Aoqpt sodel in 

which dl/dt is related to divided by the depth of the wetting front. 

* • 

For both processes, the effect of suction is represented by the but 
the aain difference is that for infiltration is considered fixed for any 
soil type, whereas for evaporation is dependent on e^. In (21) the 
evaporation rate is related to the averaged suction gradient, as represented 
by This relationship reflects the fact that bulk evaporation is 

affected by the noisture distribution in the whole zone of drying, from the 
surface down to z^. This formulation contradicts the conventional wisdom that 
says that the evaporation rate is limited by a thin, dry, highly resistant 
layer at the surface. 

The mathematical manipulations using Campbell’s parameters which yielded 
(20) also yield an explicit expression for the W factor: 

W = 6/1(l+4)tt], 

This weighting factor reflects the fundamental fact that desorption is inher> 
ently slower than sorption. Gardner (1959) arrived at the same conclusion 
based on his coxq>arison of the moan weighted diffusivity for the two processes. 

Future research may show that the concepts behind the Green-Ampt model may 
actually be more appropriate to evaporation than to infiltration. Of course, 
such a statement is mostly speculative; but for infiltration there are a 
variety of problems associated with the fact -- or at least the assumption — 
that the soil system is saturated. For instance, the effects of entrapped air 
lead to enpirical adjustments to and in the Green-Asqpt model. For unsat- 
urated systems, no such adjustments would be necessary. Nevertheless, there 
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are problems specifically associated with evaporation, for instance, the effect 
of vapor diffusion. This problem is discussed following the model test 
reported next. 



CHAPTER 3 


MODEL VALIDATION 

The CSM was tested using the same data used to develop the desorptivity 
model in Paper I. Those data came from a series of experiments conducted 
by the staff of the USDA Agricultural Research Service near Phoenix, Arizona; 
and they are described by Jackson (1973) and Jackson et al. (1976). In each 
experiment a test plot and a nearby lysimeter were intially irrigated with 
about 10 cm of water; thereafter evaporation was measured every half hour. 

The experiments were performed at different times of the year so that 
essentially only climatic conditions varied among the experiments. For the 
March experiment only, the moisture profile was measured regularly from samples 
gathered in the test plot. In testing the CSM, the average PE rate and the 
parameters of the redistribution function were evaluated directly from the 
data. Although this inforroaiton cannot be known a priori, its use allows a 
more specific test of the model. Moreover, predictive equations for these 
variables are available. 

There is one important difference in the information used to develop 
the desorptivity model and that used to test the CSM. The desorptivity model 
requires an independent estimate of the time of transition which is used 
directly in the calculations. For the CSM, this variable was used only to 
determine the average PE rate from the measured evaporation. In addition, 
the PE rate for the December experiment used in this test was 2.6 mm day 
instead of 2.1 mm day”^ as listed in Paper I. The rate was revised becatise 
subsequent analysis of the measured evaporation, reported later, suggested 
that the transition to stage II occurred on the sixth day following irrigation, 
rather than on day 10 as reported by Idso et al. (1974). 
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3.1 Results 

As shown in Figure 3 and Table 1, the CSM yielded excellent estimates of 
the bulk evaporation measured by the lysimeter. For all 4 experiments the 
largest error in E after 14 days was less than 5%. Definitive statements 
concerning the model's predictive accuracy are not justified because informa- 
tion normally not available a priori was used in the calculations. However, 
the agreement between model and measurement leads to the conclusion that after 
the transition to soil-limited evaporation, climatic factors -- wind, heat, 
temperature — are unimportant to the determination of bulk evaporation. This 
conclusion was also reached using the desorptivity model in Part I, and it 
specifically applies to stage II of evaporation. However, as discussed later, 
this conclusion must be qualified somewhat because the effect of weather on 
Stage-Ill evaporation is still uncertain. 

With respect to the time of transition, the model was fairly accurate. 

For all 4 experiments the differences between the calculated and observed 
times were less than 24 hours; this level of accuracy is the best that can 
be ejqiected since the CSM applies only to average daily conditions. Because 
diurnal variations are not simulated the model can predict a transition to 
the soil-limited stage during the night, as was the case for the experiments 
of September, March and December. Considered on an hour-by-hour basis, soil- 
water evaporation is driven mainly by solar radiation. Hence, in reality 
the transition occurs during the daytime when the available energy and the 
diurnal PE rate are greatest. 



evaporotion, E (mm) 
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The results in Table 1 show that clinate, through the PE rate« has a 
significant effect on the tine of transition. Cos^aring two different 
clinatic regimes and assuaing other factors to be identical, the climate 
with the larger PE rate yields the earlier time of transition, as one 
would e3q>ect. But wore importantly, the larger PE rate also yields the larger 
cumulative evaporation at all times, and the CSM predicts this pattern. It 
predicts this pattern even though the calculated evaporation is independent 
of PE after the transition, thus the results generated by the continuous 
similarity model lead to the same conclusions reached in Part I using the 
desorptivity model. The conclusion is that the conventional concepts of evapo- 
ration that lead to the practice of scaling actual evaporation rates to PE 
rates are unfounded from a physical basis; yet the result -- increased 
evaporation for increased PE rates -- is justified. 

3.2 Calculated Moisture Profiles 

For four selected days during the March experiment, the simulated profiles 
are plotted in Figure 4 along with soil moistures measured at noon on the 
respective days. For days 5, 9 and 14, the simulated profiles are consistently 
drier than the actual profiles. This deviation conforms to the observations 
reported by Jackson, et al. (1976). They noted that for several days following 
the transition to stage II, the surface of the test plot was invariably wetter 
than that of the lysimeter. . This is the problem of spatial variability that 
led Jackson, et al. (1976) to identify a transitional stage between the exclu- 
sively climate-limited and soil-limited stages. Because the focus of this test 
is to calculate the evaporative loss from the lysimeter, where the stage- I 
evaporation rates are precisely known, the discrepancies in Figure 4 for those 
3 days are not considered to be is^ortant to this study. 



Depth Z (cm) 
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symbols indicate data used to generate the average 
PE rate. 
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For the March experiment, additional measuremc'nts were made 15, 16, 23 
and 37 days following irrigation. For day 37, there is a large difference 
in the simulated and observed profiles in Figure 4, and this difference is 
attributable to the vapor diffusion in the near-surface region that charac- 
terizes stage III. Vapor diffusion causes increased diffusivity at low 6, 
and theoretically, it results in a "drying front" or an inflection in the 
moisture profile (van Keulen and Hillel, 1974). Although this inflection 
point is not apparent in Figure 4, clearly neither the simple exponential 
diffusivity function nor the moisture profile function (Eq. 16) required 
by the eSM is applicable under thise conditions. Consequently, it would 
be advantageous to know when the transition is likely to occur in order 
to evaluate the reliability of the model results. This problem is addressed 
in the discussion section. 

The moisture profile is also described by the simulation variables E* 
and z^. In Figure 5, the computed E* increases essentially linearly during 
stage I until the transition when the rate of increase diminishes and 
approaches zero. For all experiments, E* is approximately 1/3 of the cumu- 
lative evaporation, E, the difference being caused by the redistribution 
effect, is directly proportional to E*, and it approaches 32 cm in all 
of the experiments. 

As for the actual E* for the March experiment, it is consistently less 
than the calculated value during the f^rst two weeks, reflecting the tendency 
of the test plot to have a wetter surface than did the lysimeter. Late in 
the experiment the observed E* becomes large due to vapor losses. 




time (days) 

FIGURE 5. Model variables. The bold lines indicate E* and zj on the appropriate axes. Symbols 
indicate the observed B* for March* and the hatched areas indicate the depth interval 
of maximun observed 0 where the actual depth of drying probably occurred. 
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To illustrate the usefulness of the aoisture profile function, the 
stored soil wctcr, calculated and observed, are con^ared in Table 2. The 
stored water is calculated by integrating the noisture profile function 
(Eq. 16) between the surface and any selected depth, up to z^. In Table 2, 
the depth of integration is linited by the increments of measurements of 
the actual moisture profile. To estimate the distribution of soil water 
stored below z^, one needs a redistribution model to estimate the advance 
of the unsaturated wetting front. In Table 2, the computed results are 
reasonably accurate, but they do show the bias towards a dry profile 
relative to the field measurements during the first two weeks, followed 
by the reversed trend later due to vapor diffusion. The extrapolation of 
surface measurements to estimate soil water at deeper depths is a key problem 
in the development of remote sensing methods, and Eqs. (14) and (16) of CSM 
should prove useful. 
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TABLE 2. IfAtnr Stored BtnMth fti* Surfaet 


tlM 
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•toragt 
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W 

(») 

M 

(») 

(■>) 

% 

S.67 

10 
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25.8 

-3.2 

-12.4 

4.67 

20 

47.0 

50.9 

-3.9 
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5.67 

f< 

45.7 

48.4 

-2.7 

- 5.6 

6.67 


44.7 

46.6 

-1.9 

- 4.1 

7.67 

«i 

43.8 

45.2 

-1.4 

- 3.1 
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30 
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-1.9 

- 2.7 

9.67 

II 

66.9 

68.6 

-1.7 

- 2.5 

10.67 

•« 

65.1 

67.3 

-1.2 

- 1.9 

11.67 

If 

65.3 

66.7 

-1.4 

- 2.1 

12.67 

If 

64.6 

65.6 

-1.0 

- 1.5 

13.67 

If 

64.0 

64.4 

- .4 

- .6 

14.67 

If 

63.4 

63.5 

- .1 

- .2 

15.67 

If 

62.8 

62.8 

0.0 

0.0 

22.67 

II 

59.8 

57.1 

2.7 

4.7 

56.67 

II 

56.1 

48.9 

7.2 

14.7 


CHAPTER 4 


DISCUSSION 


At this point, the first two objectives of this report have been 
•ccoe^lished: the continuous siailerity aodel has been fully described, 

and it has been tested. The discussion below addresses the remaining two 
objectives: to exaaine the evaporation data to see if they do indeed indicate 

the stages of evaporation, and to relate observed stages of evaporation to 
the measured surface moisture content. Note that the relationships observed 
below are more tentative than the results presented so far. This is so because 
there are very few data under very dry conditions, and because data from 
the lysiaeter are related in this analysis to surface measurements made in 
the test plot, even though the previous section indicated discrepancies in 
these data sources. Nevertheless, the relationships are interesting enough 
to warrant exaainination. 

4.1 Stages of Evaporation 

To examine the bulk evaperation, it is useful to consider the log* log 
transformation of dE/dt versus E. For the a^iel results, this truisformation 
is straightforward, yielding the continuous curves shown in Figure 6. For 
the field data, the hour-by*hour variation ir the evaporation rate was 
eliminated by plotting the logarithms of each daily evaporative loss versus 
the cumilative loss at the end of that day. 

For the March experiment, data for days 23 and 37 are also included 
although there is some uncertainty in the cumulative evaporation for those 
times. The cumulative amowts were determined by differencing the moisture 
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stored in the soil colionn and by assuning that no drainage occurred. The 
estiaated range of error given this approximation is considerd to be 
negligible. 

With respect to the simulated results shown in Figure 6, the transition 
between stages T and II is clearly depicted by the abrupt change in slope; 
and the field data also exhibit this abrt^t change. In fact, the change in 
log dE/dt observed in the data between days 5 and 6 during the December 
experiment in^lies that stage II probably began sometime during day 6. This 
is the basis for adjusting the PE rate in the simulation. Additional, although 
indirect, evidence for a transition on day 6 is provided by the model itself 
which calculated the time of transition for the other three experiments. 

As for the transition to stage III, only the data for the March experiment 
exhibit a new trend under very dry conditions. I interpret the deviations 
between the simulated and observed results for days 15, 16, 23 and 37 as stage- 
III evaporation. At the end of 37 days the actual evaporation is estimated 
to be about 50 mm, whereas the calculated amount is only 37 mm. Thus, the 
CSM is not applicable for prolonged periods of stage-III evaporation. Further- 
more, the pcssible effects of climate on stage III are simply unknown. 

As for the other experiments, I infer that the slight increases in log 
dE/dt for day 5 during July and for day 8 in September may result from stage- 
III condtions. This inference is based on the measurements in the surface 
moisture content described next. In any event, the bulk evaporation rates did 
not deviate much from the simulated results so there is no clear evidence of 
stage III during the July and September experiments. 


tvaporofion rott, dC/dt (mm doy'^) 



FIGURE 6. Transformed evaporation data. Each data symbol corresponds 
to one day, and certain symbols are individually identified 
by the number of days following irrigation. Simulated re> 
suits for July, September, March, and December extend to 
7, 14, 37 and 14 days, respectively. Circled symbols 
indicate data used to generate the average PE rates. 
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4.2 Surface Soil Moisture 

Past researchers have stated that the surface Boisture content can serve 
as an indicator of the different states of evaporation. For instance, Hillel 
(1971) stated that stage II begins when the surface is approxiaately air-dry. 
In contrast, Jackson (1973) stated that stage II begins when the surface has 
an intexaediate wetness and that stage III begins when the surface is air-dry. 
With these conflicting views in aind, the piurpose of this subsection is to 
relate the observed stages of evaporation to aeasureaents of 6^. 


Continuous aeasurraents of the near- surface aoisture profile were aade 
only during the July and March experiaents. I investigated the average daily 
soil aoisture, 6^, and fmmd trends in log 6^ versus time that roughly 
correspond to the stages of evaporation identified above. For March, 

Figure 7 shows three linear trends in log 6^, with the suggestion of a fourth 
one at the beginning of the experiaent. Interestingly, the trends in log 
are alaost perfectly linear. The physical swaning for this linearity is unclear, 
but I aa presently investigating these relationships. 

In Figure 7, is approxiaately .40 ca^ ca~^ at t ■ 0, as extrapolated 
froB deeper aeasureaents. Assuaing a straight-line trend, the first breakpoint 
in the data corresponds to the transition to stage II. In keeping with the 
observation of Jackson (1963), stage II begins when the surface has an inter- 
aediate 9^. The second breakpoint occurs in the aiddle of stage II and corre- 
sponds to an average surface aoisture of about .06 ca^ ca'^ which is the air- 
dry value of 9 for the particular soil, Avondale loaa (Jackson, 1963). The 
break point probably signifies the beginning of interaittent vapor diffusion 
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in the near-surface region, but there is no corresponding change in the 
evaporation rate. Thus, bulk evaporation is still United by the novement of 
liquid water within the profile, and I have labeled this tine interval as 
substage Ilb. A third breakpoint occurs about day 14 when stage III began. 

In fact, day 14 is the first day when 6^ remained less than air-dry through- 
out the day. Consequently, I conclude that stage III can connence when the 
surface remains less than air-dry throughout the entire day, and perhaps this 
"sub air-dry" condition could serve as a criterion for stage III. 

There are few data for July, but based on the information for March, the 
saro system of trends can be imposed on the data that are available. The 
data indicate that the surface became sub air-dry on day 5, much sooner than 
in March. However, there is no evidence of stage III in the evaporation data, 
so it appears that the soil surface may become sub air-dry while the evaporation 
rate declines in stage-II fashion. Therefore, a sub air-dry surface is a 
necessary, but not sufficient, condition for stage III; and it cannot be 
used as a sole criterion for stage III. 


CHAPTER 5 


CONCLUSIONS 

Idike et tl. (1977) deaonstrated that infiltration under field 
conditions can be calculated fron the ainioal parameters required to 
describe the soil's hydraulic prq[>erties. Analogous to the Green-Aaq>t 
infiltration equation, the continuous siadlarity model describes evapo- 
ration from an unvegetated soil; and it, too, depends on fundamental 
hydraulic parameters. As doaonstrated herein, the C^ is applicable under 
field conditions. 

The fact that the model performed well in all 4 of the experiments 
tends to support the assuaqptions iq;>on which the model is based. Based on 
the success of the simulations, it is concluded that for an initially moist 
soil, redistribution can be calculated independently of evaporation. Indeed, 
a submodel for redistribution must be incorporated if the CSM is to be a 
predictive model. 

It is also concli^ed that bulk evaporation during stage II can be 
calculated without regard to the surface soil moisture. 6 q does not affect 
dE/dt because of the nonlinearity of the fundamental hydraulic properties of 
soil. When is moderately low the moisture gradient beneath the surface is 
steep. If 6 q is decreased, the diffiisivity decreases, but the gradient 
increases so that dE/dt remains virtually unchanged. Stated in a different 
manner, the soil water pores represented by the low range of 6 never dry to 
any appreciable depth during stage II. 
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However, this conceptualization is not true during stage III, and it 
mist be slightly Bodified when diurnal variations are considered, as was 
shown by Jackson (1963). With the daily pulses of sunlight, the surface 
layer also acts like a tenporary storage coiqiartnent , discharging aoisture 
during the day that accumilated fraa deeper with the profile during the night. 
However, the results reported here indicate that this storage effect does not 
control bulk evaporation. The controlling factor is the 6 gradient between 
the surface and the depth of drying. In the nodel this gradient is represented 
by the state variables 6^ and E*, where E* is directly related to the depth of 
drying. 

This conceptualization differs with the often-stated viewpoint that 
evaporation is limited by a thin, dry, highly resistant layer at the surface. 
In fact, when the surface becomes air dry, its diffusive capability increases; 
and it can act as a wick, offering relatively little resistance to moisture 
moving upward from within the soil. 



APPENDIX 


COMPUTER CODE FOR THE CONTINUOUS SIMILARITY MODEL 
OF SOIL-WATER EVAPORATION 


The FORTRAN prograa is listed on the next two pages^ followed by a listing 
of the program's output. The program contains all the necessary information 
in the DATA statements, so input is not required. The parameters of the 
functions for diffusivity and redistribution are output using the NAMELIST 
convention. Subsequently, the simulation is con^leted four times, once for 
each of the seasonal e]q>eriments. The solution is obtained using a second- 
order, Runge Kutta integration method. The method requires that the rate 
equations for the evaporation deficit and cumulative evaporation (E* and E, 
respectively) be evaluated twice at each time step. The rate equations are 
computed in the subroutine RATE which conforms to the flow chart given in 
Figure 2. The output, including values for the state variables and for their 
instantaneous rates, is printed for each day of the simulation. 
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